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Abstract 

We derive exact strong-contrast expansions for the effective dielectric tensor e e of electromag- 
netic waves propagating in a two-phase composite random medium with isotropic components 
explicitly in terms of certain integrals over the n-point correlation functions of the medium. Our 
focus is the long- wavelength regime, i.e., when the wavelength is much larger than the scale of 
inhomogeneities in the medium. Lower-order truncations of these expansions lead to approxima- 
tions for the effective dielectric constant that depend upon whether the medium is below or above 
the percolation threshold. In particular, we apply two- and three-point approximations for e e to a 
variety of different three-dimensional model microstructures, including dispersions of hard spheres, 
hard oriented spheroids and fully penetrable spheres as well as Debye random media, the random 
checkerboard, and power-law-correlated materials. We demonstrate the importance of employing 
n-point correlation functions of order higher than two for high dielectric-phase-contrast ratio. We 
show that disorder in the microstructure results in an imaginary component of the effective dielec- 
tric tensor that is directly related to the coarseness of the composite, i.e., local volume-fraction 
fluctuations for infinitely large windows. The source of this imaginary component is the attenua- 
tion of the coherent homogenized wave due to scattering. We also remark on whether there is such 
attenuation in the case of a two-phase medium with a quasiperiodic structure. 
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I. INTRODUCTION 



The problem of determining the effective dielectric tensor and other mathematically 
equivalent properties of random media dates back to the classic work of Maxwell [1]. Cal- 
culation of the effective dielectric tensor of disordered composite materials is essential for a 
wide range of applications, including remote sensing (e.g., of terrain, vegetation, etc.) 
the study o 



composites 



' wave propagation through turbulent atmospheres [3|, active manipulation of 
4(, and as a probe of artificial materials |5j. The effective dielectric tensor at 
ong wavelengths plays a particularly important role in the study of electrostatic resonances 

This paper is concerned with the calculation of the effective dielectric tensor e e of a 
two-phase dielectric random medium associated with electromagnetic wave propagation in 
the long wavelength regime, i.e., when the wavelength is much larger than the scale of 
inhomogeneities in the medium. The complementary regime, in which the wavelength is 
much smaller than the inhomogeneities, may be studied numerically using a ray-tracing 
technique based on geometrical optics The purely static (as opposed to the full dynamic) 
problem, on which Maxwell's work was centered, can be considered to be a special case of the 
present work in the limit of infinite wavelength (or zero frequency) . Although there has been 
extensive previous work on the dynamic problem (see Refs. andjic] and references therein), 
the vast majority of studies that attempt to relate e e to the microstructure emp 



point correlation information 



ion 



121 ] (with a few exceptions, see, e.g., Ref. 



oy 



only two- 



13). On the 



other hand, there has been a great deal of work on the static problem incor por atin g th ree 
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point and higher-order correlation functions; see, for example Refs 
and references therein. 

Here we present, for the first time, explicit closed-form series expansions for the effective 
dielectric tensor of two-phase random media in three dimensions in terms of certain integrals 
of n-point correlation functions for the dynamic problem in the long-wavelength regime. 



The approach used follows directly from one given origina 



19 



ly for the purely static problem in 
. An advantageous feature of this 



any dimension in Ref. [17J, and expanded upon in Ref. 
formalism is that it gives rise to expansion parameters involving the dielectric constants that 
yield very good convergence properties even for high phase-contrast ratio. The technique 



is called the strong- contrast expansion for this reason. Tsang and Kong [IJj have employed 
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a similar method, but only up to the two-point level and they provide no justification for 
or limitations on the class of microstructures and phase-contrast ratios for which their two- 
point approximation is applicable. We restrict ourselves to considering component materials 
that are themselves isotropic but that generally possess complex-valued dielectric constants. 
As we will demonstrate, lower-order truncated forms of these expansions can serve as useful 
approximations of the effective dielectric tensor. Although the formalism is applied here 
to electromagnetism in three dimensions, it is straightforwardly generalizable to a class of 
other vector fields in arbitrary dimension 19]. 

We obtain new approximations for e e by truncating the exact strong-contrast expansions 
at the two-point and three-point levels that depend upon whether the medium is below or 
above its percolation threshold. Using these two-point approximations, we estimate e e of 
dispersions of hard spheres, hard oriented spheroids and fully penetrable spheres as well as 
Debye random media, random checkerboard, and power-law-correlated materials. We use 
a three-point approximation to evaluate the effective dielectric constant in the case of the 
fully-penetrable-sphere model in order to establish the importance of three-point information 
at high phase-contrast ratio and volume fraction, and assess the accuracy of the two-point 
approximations. For many of the model microstructure considered here, the explicit forms 
of the corresponding rz-point correlation functions follow from the general representation 
formalism of Torquato and Stell (2^] . 

A significant qualitative difference between the purely static problem and the dynamic 
problem is that in the latter, the effective dielectric tensor may be complex even if the 
component materials have purely real dielectric constants. This results from incoherent 
scattering of the incident wave. The imaginary component of the effective dielectric tensor 
has special significance to remote sensing applications, as it plays a key role in, for example, 
the calculation of backscattering coefficients 11]. Although the expansions we derive here 
are applicable in general to materials with complex dielectric constants, we will apply them 
to cases in which the material phases have real dielectric constants. We do this because 
in this regime, disordered media will in general yield nonnegative imaginary contributions 
to the effective dielectric tensor purely due to scattering, not absorption. The physical 



manifestation o 
experiments 21 



this can be understood in the context of optical or ultrasonic transmission 



22] , in which a wave pulse propagates through a finite-sized slab composed 



of a disordered material, and the transmitted signal is measured. The transmitted field 
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is observed to be a superposition of a coherent pulse representing the incident signal and 
uncorrelated noise due to random scattering. The attenuation of the incident pulse is a 
consequence of this scattering. Since periodic media propagate waves without any loss (and 
thus with purely real effective dielectric tensors), they do not exhibit imaginary parts in their 
effective dielectric tensors. In the case of statistically homogeneous and isotropic media in 
particular, we show here that the leading-order contribution to the imaginary component of 
the effective dielectric constant is directly related to the coarseness of the composite for 
an infinitely large window. The coarseness is defined as the standard deviation of the volume 
fraction of a given phase, in some observation "window" within the composite divided by the 



total volume fraction of that phase [19|, [23|. It has recently been sugges ted that local density 



fluctuations provides a crude measure of disorder of a system |24j . [25] . In the language of 
Ref. |24|, a hyperuniform system is one in which = 0, which has been shown to always be 
the case for periodic systems, and therefore is consistent with the fact that Im [e e ] = for 
such media. Note that a non-zero imaginary component of the effective dielectric constant 
can be directly associated with a mean free path of waves propagating through the given 
medium 10]. 

In Section II, we introduce our formalism, derive the strong-contrast expansion for the 
effective dielectric tensor, and discuss the approach towards applying these techniques at 
the three-point level. In Section III, we discuss a variety of model microstructures not pre- 
viously studied in the present context and their associated two-point correlation functions. 
In Section IV, we present results for the model microstructures discussed in Section III, 
showing volume fraction and dielectric-contrast ratio dependence at the two-point level, and 
demonstrating the importance of three-point calculations at sufficiently high phase contrasts 
and volume fractions. In Section V, we discuss conclusions based on our results. Appendix 
A presents a generalization of our formalism, up to the two-point level, wherein the expan- 
sion is carried out with an arbitrary "reference" or "comparison" material (as described in 
Section II. A), discussed hence). Appendix B gives a short proof of why periodic media, at 
low frequencies, must give rise to effective dielectric tensors with no imaginary component. 
Finally, in Appendix C, we simplify the key two-point integral for statistically homogeneous 
but anisotropic media with azimuthal symmetry, an example of which is a dispersion of hard 
oriented spheroids. 
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II. THEORY 



Here we extend the formalism for determining strong-contrast expansions of the purely 



13, 



19| to the dynamic case in 



static effective dielectric tensor in arbitrary space dimensions 
M 3 when the wavelength is much larger than the inhomogeneity length scale. This method 
is based on finding solutions of certain integral equations using the method of Green's 
functions. 



A. Strong- Contrast Expansions 

Following Ref. H9 



(which gives greater detail than Ref. [17|), we begin by considering a 
macroscopically large ellipsoidal sample composite material in M 3 composed of two phases, 
labeled 1 and 2, which is itself embedded in a homogeneous "reference" (or "comparison) 
material of dielectric constant Eq. We choose a specific macroscopic shape to call attention 
to the fact that the average fields in the problem are shape dependent. An advantage of this 
formalism is that it eliminates the shape dependence of the effective dielectric tensor. The 
length scale of inhomogeneities within the composite is assumed to be much smaller than 
the shape itself. The two phases have dielectric constants E\ and 62, and we define indicator 
functions in the following way: 

, s I 1, if r lies in phase p 

2>)(r) = I (1) 
I 0, otherwise 

for p = 1,2. Thus, the volume fractions of the phases are (Z^(r)) = <p p and <p q = 1 — <p p 
(p / g). We may therefore write e(r) = SiX^{y) + e 2 X( 2 )(r) i n the composite. Note that 
the dielectric constants S\ and Ei may be complex. 

When solved for the electric field, time-harmonic Maxwell's equations reduce to 

VxVxE(r)- £ (r)(-) 2 E(r)=0, (2) 



c . 

with E the electric field, six) the reduced dielectric constant, ui the frequency of the time- 
harmonic solution, and c the speed of light. We assume here that \xj ' \i§ = 1, and that the 
component dielectric materials are themselves isotropic. We can rewrite this homogeneous, 
linear equation in a form that is suggestive of perturbation theory: 



V x V x E(r) - M "J E ( r ) = ( £ ( r ) " e ^ {-) E ( r )> ( 3 ) 



where we have taken the comparison material to be one of the phase materials, i.e., Eq = e q , 
where q is either 1 or 2. We employ this choice for simplicity here, but this is not essential; 
in fact, depending on the details of the structure in question, it may improve convergence to 
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271 ] . We discuss this point in further detail 



choose a different comparison material 
in Appendix A. 

As a shorthand, we will define kg = o q = (^) 2 e q , for any phase, q = 1 or 2. In order to 
solve this equation for an arbitrary structure via perturbation theory, we require a Green's 
function for the operator on the left side of Eq. (j3J), which must therefore be a tensor. This 
is the dyadic Green's function. In three-dimensional spherical coordinates, it is given by 

G(r, r') = -^-5(r - r') + Gi(r - r')I + G 2 (r - r')rr, (4) 
da q 



where 



d(r - r') = (-1 + ik q r + k q r 
G 2 ( r — r ') = (3 — 3ik q r — a q r 



47r/c^r 3 ' 

2 )~ 

Arck^r 3 



with r = |r — r'|, r is a unit vector directed from r' towards r, and I is the unit tensor in 
three dimensions. In the limit k q — ► 0, we recover the static result [3] in three dimensions. 
Note the delta function in this expression; this is the dipole "source" of the radiation; its 
coefficient is dependent on the shape of the "exclusion volume" around this source. For the 
coefficient shown, the exclusion volume must be spherical in shape. 

The Green's function satisfies the following partial differential equation: 

V x V x G(r, r) - a q G(r, r') = 16 (t - r'). (5) 

This implies that we can write the electric field as 

E(r) = E (r) + J dr'G(r, r') [a{v) - a q \ E(r'). (6) 

We can express this integral equation more compactly in linear operator form: 

E = E + GP, (7) 

where we define the polarization vector field as 

P = [a(r) - a q ] E. (8) 



The next step is to extract the delta function contribution from the Green's function solution. 
In doing so, we obtain a new field, F, the cavity intensity field, which is directly proportional 
to E. The resulting integral equation is 

F = E + HP, (9) 

where H is the principle value of the Green's function in Eq. (j3J), namely 

H(r - r) = Gi(r - r')I + G 2 (r - r')rr. (10) 

Here, we have defined 



I 4 



da q 



F(r) 

where I is the unit dyadic tensor. Now, P and F are directly related: 

air) - a, 



E(r), (11) 



P(r) = jM- i -F(r). (12) 

or, implicitly defining L, 



1 + ^< r > 



P(r) = L(r) ■ F(r). (13) 

Note that L is an isotropic tensor. Instead of using the standard definition of the effective 
dielectric tensor, (D) = e e - (E), we may equally well use the above equation. Thus, implicitly 
defining L e , we have 

<P(r)) = L e -(F(r)>, (14) 

where (.) denotes the ensemble (volume) average (assuming ergodicity). We may write L e 
explicitly at this point. It is given by 



a P — o ,, 



L e = L e I = ^ q - , (15) 

where a e = (^) 2 £ e , with e e the effective dielectric tensor of the composite. 

We now have an equation that defines the effective dielectric tensor. The last step involves 
eliminating the background field, Eo. This is done because, as is known from electrostatics, 
the relationship between the applied field and the average fields (P and L, for example) 
is shape-dependent. Thus, eliminating E results in an effective dielectric constant that is 
independent of the shape of the macroscopic ellipsoid. It has been demonstrated that the 
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elimination Eq thus leads to integrals of correlation functions within the formalism that are 



absolutely convergent 



17 



19] . 



We define here a new tensor operator S in the following way: 

S = Lfl-LHl \ 



(16) 



It follows directly that (P) = (S)E , and we may thus eliminate E from Eq. (J9]). Upon 
taking the ensemble average of the latter, we may write (F) in terms of (P), the first few 
terms of which are given by: 



(F(rO> 



dr 2 dr 3 



(£(n)) 



(L(ri)L(r 2 )) — (I/(r 1 ))(L(r 2 )) 



<L(n)><L(r 2 )> 



H(r lf r 2 )(P(r 2 )>- (17) 



(L( ri )L(r 2 )L(r 3 )) (L( ri )L(r 2 ))(L(r 2 )L(r 3 )) 
(L( ri ))(L(r 2 )) (L( ri ))(L(r 2 ))(L(r 3 )) 



H( ri ,r 2 )H(r 2 ,r3)(P(r 3 ))-.. 



(18) 

From this expression and Eq. (JHJ), we obtain the following exact expansions for statistically 
homogeneous media involving the effective dielectric tensor: 



(5 2 pq f p {a e - a q l)-\a e + 2a q l) = 4> p (l pq l - ^ A<*%, (p ^ q), 



where p ^ q, 



n=2 



6p Eq 



pq 



(19) 



(20) 



a p + (d— l)a q e p + (d—l)eq 
and the n-point tensor coefficients A„ are certain integrals over the n-point correlation 
functions S^f 1 associated with phase p. In particular, for n = 2 



A 



(p) 



d_ 



dr&H 



where 



t (p) (r) = nk 2 „H(r) 



(21) 



(22) 



and Q is the solid angle of a sphere in dimension d, and H(r) is given by Eq. (flQl) . For any 
n > 2 , we have 



y <Zr 2 ... / rfr n t^(r 1; r 2 )-t«(r 2 , r^-t^r^, r n ) A^(n, r n ), 

(23) 



9 



where A^(ri, ...,r n ) the determinant is given by 



A^ 



4 p) ( ri ,r 2) r 3 ) Si p \r 2 ,r 3 ) 








Sn\ri, r n ) S , ^ 1 (r 2 , r n ) ••• S^^r^i, r 



7 (p) 



(*>), 



(24) 



and 



5«(n, r 2 , . . . , r„) = (/W( ri )/W(r 2 ) • • • J«(r B )>. (25) 

is the n-point correlation function of phase p. For statistically homogeneous media (the case 
for which this formulation applies), the latter quantity is translationally invariant and there- 
fore depends only on relative displacements, i.e., Sn ( r i, r 2, • • • , r «) = 5n (i"i2> r i3, ■ ■ ■ , r in)i 
where we have chosen point 1 to be the origin, and, in particular, (ri) = P . 
Remarks 

1. We see that the integral-equation approach given here is entirely general. It is equally 
well suited to any problem in which the local and averaged constitutive relations have the 
same form, and thus the appropriate dyadic Green's function would replace Eq. (J4]). Indeed, 
the determinant fl24|) for the dynamic problem considered here is exactly the same as that 
in the static problem, as given originally in Ref. Il7. 



2. Note that Eq. (fT9|) represents two different series expansions: one for q = 1 and p = 2 
and the other for q = 2 and p — 1. 

3. The quantity j3 pg , given in Eq. fl20|) . is the strong-contrast expansion parameter, the 
form of which is a direct consequence of the choice of exclusion volume associated with the 
source-term of the Green's function. Any shape besides the sphere would have necessarily 
led to a different expansion parameter and therefore to significantly different convergence 
properties jl9|] . Clearly, (3 pq may lie within the range — (d — 1) _1 < (3 pq < 1. The radius 
of convergence of the series given in Eq. (|T9l) is thus greatly widened beyond that of a 
weak-contrast expansion (i.e., the simple difference of the phase dielectric constants). 

4. In the purely static case, the series represented by f|T9|) may be regarded as expansions 
that perturb around the optimal structures that realize the generalized Hashin-Shtrikman 



bounds 



derived by Willis [29(], as discussed by Torquato [19]. In particular, these optimal 



structures are certain dispersions in which there is a disconnected, dispersed phase in a 
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connected matrix phase. The lower bound corresponds to the case when the high-dielectric- 
constant phase is the dispersed, disconnected phase and the upper bound corresponds to 
the instance in which the high-dielectric constant phase is the connected matrix. Thus, we 
expect that for the dynamic problem under consideration, the first few terms of the expansion 
(TE?j) with q = 1 and p = 2 will yield a reasonable approximation of e e for two-phase media, 
depending on whether the high-dielectric phase is below or above its percolation threshold, 
as discussed further in Section II. C. 



B. Macroscopically Isotropic Media 



Consider, as we will to a large extent in this paper, two-phase media that are statistically 
homogeneous and isotropic, and thus macroscopically isotropic. In this case, the effective 
dielectric tensor is proportional to the identity tensor, and can thus be treated as a simple 
scalar, namely, the effective dielectric constant of the medium. This of course simplifies the 
calculation enormously. We thus take the trace of both sides of Eq. ( 1T91) and divide by d. 
Thus, the full expression in the isotropic case is given by 



where A 



(p) 



Tr 



(p) 



n=2 

I d and (3 eq is the effective polarizability, given by 



(26) 



eq 



e p +(d- l)e q 

The scalar two-point coefficient, as specified by taking the trace of (J2T1) . is given by 



(27) 



.4 



(p) 



d 

n 



dr 



Tr [t 



(p)i 



d 



S. 



(p), 



r - 



2k 2 / drexp(ik q r)r 



S. 



(p), 



r — 



(28) 



with this integral being straightforwardly carried out either numerically or analytically, 
depending on the form of the correlation function. Provided that the correlation function 
S% (r) decay sufficiently rapidly to its long-range value of the integral in (!28~i) will be 
convergent. 

Note that by Eq. ( J28l) . A% must be zero in the purely static problem (a; = 0) if the 
medium is statistically homogeneous and isotropic. Thus, in the static problem, assuming 
statistical homogeneity and isotropy, two-point information is actually incorporated even 
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This suggests that our 



though A2 is zero. This subtle point is elaborated in Ref. [19 
formalism is extremely well suited to the non-static problem in the low frequency limit 
being considered here. Expanding the two-point coefficient A% given in Eq. fl28l) through 



third order in k q about k q = yields 



A? 



2ki 



drr 



S. 



(p). 



+ 2ikl 



drr 



(p), 



+ 0(kt). (29) 



It should be expected that the imaginary component of the effective dielectric constant of 
statistically homogeneous and isotropic systems should be positive; otherwise the homoge- 
nized coherent wave would be amplified rather than attenuated. A simple analysis of Eq. 
(|29|) bears this out. The second term on the right hand side, i.e., the leading-order imaginary 
component, is nothing more than the zero- wave- vector structure factor of the composite ma- 
terial. Since s!f\r) is obtained from a realizable configuration, this is necessarily positive 
30( | . This is true because the structure factor is nothing but the squared norm of the Fourier 
transform of the function X^ p \r) — (fi p . The leading-order term of the imaginary component 
in the expansion given in Eq. ( 1291) is directly proportional to the volume integral of the 



function S. 



(p), 



This is exactly (02Coo) 2 t>o, where is the coarseness of the compos- 



ite structure in the limit of an infinitely large window, and Vq is the window volume 



19 



23]. 



Thus, the value of the nonnegative imaginary part of the effective dielectric constant to 
leading order is determined by local-volume-fraction fluctuation over v ery large windows, 
which may be taken to be a crude measure of disorder in the system. [24J Note also that 
Cqo is proportional to the single-scattering intensity of scattered radiation from the random 
medium in the infinite- wavelength limit 19|, [31 1. 

Another important fact that emerges from this formalism is that at the two-point level 
the correction to the static effective dielectric constant only enters at second order in k q (i.e., 
second order in the frequency). It is hence a relatively small contribution. The imaginary 
term enters only at third order in the wave number. 

Now, let us consider the simplification of the three-point coefficient A^\ From the general 
expression ([23]) . we find 



.4 



(p) 



9 



d 3 r 



12 



rf 3 r 13 t^(r 12 )-tW(r 12 ) 



S3^( r 12> r 13) — 



(p), 



(30) 
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The first few terms of the Taylor expansion of (1301) about k q = is given by 



A 



(p) 



18 fd 3 r 12 fd 3 r 13 



(47T) 



9 



00 

r 



12 

s 



13 



-i 



6 



l + 4(r 2 + s 2 ) 



S3' ( r 12, r 13) — 



5. 



5 3 W (r, a ,Ai) 



+ 0{k\) 
i d(^). (31) 



where /1 = cos6*, 6 1 is the angle between ri 2 and ri 3 and P 2 is the second-order Legendre poly- 
nomial. The simplified second line is obtained by exploiting the homogeneity and isotropy of 
the medium, reducing the original six-dimensional integral to a three-dimensional one 18]. 
Here r = |r 12 | and s — |r 13 | are the side lengths of a triangle and 9 is the angle between these 
sides. Note that the third-order term in k q is exactly zero, the next-lowest-order contribu- 
tion to the real component must be of the order of k q (or higher), and the next-lowest-order 
contribution to the imaginary component must be of the order of k q (or higher). 

In Section IV of this paper, estimates for the effective dielectric constant that include 
three-point information will be presented for the fully-penetrable- sphere model. 



C. Two- and Three- Point Approximations 

Practically speaking, it is difficult to ascertain four-point and higher-order correlation 
functions which therefore prohibits an exact evaluation of the expansion for the effective 
dielectric tensor in Eq. (I19p . However, as shown in the static problem [171. [19|, lower-order 
truncations of the exact expansion (|T9|) at the two-point and three-point levels have proved 
to be accurate approximations for the effective dielectric tensor. For the macroscopically 
isotropic case, the two-point and three-point approximations used to calculate the effective 
dielectric constant, obtained by truncating the series in Eq. ( 1261) . are given by 



Ep Eq 



E p + (d- l)E q 



_E e + (d- l)e Q 



Ep Eg 



_E p + (d— V)E q 



ip^q), (32) 



and 



E p + (d- l)e q _ 



E e + (d- l)E q _ 



Ep Eq 



e p + (d- l)e q _ 



-A 



(p) 



e p + (d- l)e q 

(33) 

respectively, where A^ is given by Eq. (1291) . and A^ is given by Eq. ( 13TT) . The two-point 
approximation (|32|) is exact to second order in e p — e q for any <p p and is exact to first order in 



, (p ^ q), 
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(j) p for any phase-contrast ratio e p /e q . The three-point approximation (|33|) is exact to third 
order in e p — e q for any (p p and is exact to first order in <j) p for any phase-contrast ratio e p /e q . 
However, the most difficult cases to treat theoretically are when both the volume fraction and 
phase-contrast ratio are significantly different from zero and unity, respectively. The strong- 
contrast approximations (13"2"j) and (I33p are expected to provide reasonable estimates of e e 
for a certain class of dispersions in this more difficult regime because they are perturbations 
of the strong-contrast expansions in the infinite- wavelength limit (pure static case), which 
have been shown to be in excellent agreement with both precise computer-simulation and 
experimental data for a variety of dispersions 0, Q]- Specifically, if e p > e q , ( 1321) and ( |33l) 



with q = 1 and p = 2 will yield good estimates of e e provided that phase 2 is below its 
percolation threshold and that the typical cluster size of phase 2 is sufficiently small [l?], Q] ■ 



On the other hand, if phase 2 is above its percolation threshold, fl32|) and fl33|) with q = 2 
and p = 1 should provide good estimates of e e . 

III. MODEL MICROSTRUCTURES 

We study a number of model microstructures and their corresponding two-point correla- 
tion functions. In particular, we examine dispersions of hard spheres, hard oriented spheroids 
and fully penetrable spheres as well as Debye random media, the random checkerboard, and 
power-law-correlated materials. Except for the hard-ellipsoid model, all of the other models 
are statistically homogeneous and isotropic. We take a to be a characteristic length scale 
for each model. These microstructures, depicted in Fig. [1] are explicitly described in this 
section, and we take q = 1 and p = 2. The two-point correlation functions for each of these 
models is depicted in Fig. at volume fraction 2 = 0.5, with correlation lengths roughly 
equal to one another. The "correlation length" l c is the range [0, l c ] in r over which the 



magnitude of s!p — <p\ is appreciably fluctuating around zero [19]. 



A. Equilibrium Hard Spheres 

Here we consider the well-known equilibrium distribution of hard spheres of radius a 
(phase 2) in a matrix (phase 1) [34 1. For this model, depicted in Fig. all non- 

overlapping configurations are equally probable. This equilibrium model is athermal, i.e., 
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(A) 



(B) 



(C) 




FIG. 1: (Color online) Two-dimensional slices of the three-dimensional microstructures described 
in this section at volume fraction fa = 0.5. They are hard spheres (A), hard oriented spheroids 



321 ] (D), random checkerboard (E) (in 



32] (F). Both 



(B), fully penetrable spheres (C), Debye random media 

any plane perpendicular to a principal axis), and power-law-correlated materials 
the fully-penetrable-sphere model and random checkerboard at <j>2 = 0.5 are above their respective 
percolation thresholds for the black phase 2, even though planar cuts through these samples do not 
reveal that the black phase is indeed percolating in three dimensions. We expect that the Debye 
random medium and power-law-correlated materials shown here also percolate at <p2 = 0.5. 



the behavior of the system is temperature- independent. The pair correlation function g2{r) 
for particle centers may be expressed in closed form in Fourier space using the Percus-Yevick 

(2) 

approximation, as described in Ref. [19|. The two-point function S 2 (r) can be obtained from 



the pair correlation function via the following equation 



35j: 



S^\r) = pv2 %nt, {r; a) + p 2 g 2 <E> m <8> m, 



(int) 



(34) 
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FIG. 2: (Color online) Plots of the two-point correlation function Svf (r) for the five isotropic 
models detailed in this section: hard spheres, fully penetrable spheres, Debye random medium, 
random checkerboard, and a power-law-correlated material. The volume fraction for each model 
is <p2 = 0.5. Note that the correlation lengths for each model are roughly equal to one another. 



where p is the number density of spheres, Uj \r; a) is the intersection volume between two 
spheres of radius a that £1X6 cL distance r apart, which is given by 

vt t \r-,a)={ 1 ~^ + ^^ )3 ' ifr<2a , (35) 
I 0, otherwise 

the symbol <E> denotes the convolution of two functions, and the step function m is defined 
by 

I 1, if |rl < a 

m(r) = I . (36) 

I 0, otherwise 

Note that the volume fraction of spheres is given by 02 = p47ra 3 /3. The two-point correlation 
function S 2 (r) is plotted in Fig. [2j The most convenient method of obtaining Sv; (r) via 
Eq. ()34p is via Fourier transform techniques. The reason for this is that in Fourier space, 
the convolution operators become simple products, and the Percus-Yevick pair correlation 
function may be expressed analytically. In order to obtain (r) in real space, numerical 
inversion from Fourier representation must be performed. 

Of all the isotropic disordered models considered here, the hard-sphere system may be 
thought to possess the greatest degree of order because its coarseness is minimized among 
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the structures considered. For example, consider volume fraction 02 = 0.5, at which spheres 
are slightly above their freezing volume fraction. The correlation function S^\r) — 0?, 
exhibits oscillation above and below zero, suggesting strong short-range correlations and 
anti-correlation among the spheres in the system. Calculations show that as a result of this 
property, the imaginary component of the dielectric constant is very small compared with 
that of other model microstructures considered here. 



B. Equilibrium Hard Spheroids 

We also consider an equilibrium dispersion of hard spheroids (phase 2), or ellipsoids 
of revolution, in a matrix (phase 1), which are constrained to have the same orientation 
(see Fig. [D(B)) along the z axis. Because there is an axis of symmetry, spheroids possess 
azimuthal symmetry. This statistically homogeneous but anisotropic dispersion is considered 
here in order to demonstrate the application of our formalism to macroscopically anisotropic 
media, i.e., materials with an effective dielectric tensor with non-equal diagonal terms in the 
principal axes frame. The spheroid shape is defined by 

(x 2 + y 2 )/a 2 + z 2 /b 2 = 1, (37) 

where a and b are the semi-axes of the spheroid with b being along the axis of symmetry 
(i.e., z axis). We define the aspect ratio as b/a so that b/a > 1 and b/a < 1 corresponds 



to prolate and oblate spheroids, respectively. In a previous work, Torquato and Lado [36] 
showed that the correlation function for any dispersion of oriented spheroids (in equilibrium 
or not) at number density p could be transformed directly from the hard-sphere correlation 
function at the same number density by an affine (linear) transformation of the coordinate 
system. Taking phase 2 to be the spheroid phase, this transform is defined by 

Si% s (v; b/a) = S (2) HS [* o (r u /*(0)); M > ( 38 ) 

where 

<*{0) = T72> ( 39 ) 

[1 - (l-a 2 /6 2 )cos 2 ^i 1/2 



(21 

where o"o = 2a, 9 is the polar angle between the z axis and the radial vector r, and Sr, ^ 5 (r; 1) 
is the hard-sphere correlation function at the same number density. Here we employ the 
equilibrium hard-sphere model to get the corresponding expression for the equilibrium hard- 
oriented-spheroid system. 
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C. Fully Penetrable Spheres 



A fully-penetrable-sphere model is composed of spheres of radius a (phase 2) whose 
centers are completely uncorrelated in space [l9] (see Fig. [D(C)). For this dispersion, the 
particle phase percolates at 2 = 0.2895 ± 0.0005 {3^ (i.e., the medium contains infinite 
clusters of phase 2), and the matrix phase (phase 1) percolates until a volume fraction of 
spheres of 97%. The n-point correlation function for this system, for phase 1, is given in 
Ref . [l9j as 

S£\r n )=ex V [-pv n (v n ;a)], (40) 

where the number density p is defined here by <pi = exp(— r/), with rj = p4ira 3 /3, and the 
function v n (r n ;a) gives the union volume of n spheres of radius a with centers defined by 
the position coordinates r n = ri, r 2 , . . . , r n . Thus, we may write, for the particle phase, 



S?\r) 



1 - 202 + exp 



-7f- 



V2{r;a) 
vAa) 



where we have 



Vila) 



20(r - 2a) + 



3 r 1 / r \ ! 
+ 4a ~ 16 \a 



B(2a-r) 



(41) 



(42) 



and 6(x) is the unit step function. The correlation function Sr, '(r) is plotted in Fig. [21 
Here, we have used the fact that we can relate the union and intersection volume of two 
spheres by the equation 

v 2 (r;a) = 2v 1 (a)-vt t) (r ] a). (43) 

The function may be obtained in similar way. It involves an expression for the intersec- 
tion volume of three spheres, given originally by Powell in Ref. [38|, and discussed in detail 
in Ref. [l9j. In the next section, we present estimates of the effective dielectric constant 
as obtained from two- and three-point approximations for this model. Compared with the 
hard-sphere model, the fully-penetrable-sphere dispersion is significantly more disordered. 
This is reflected in the appreciably large difference between the imaginary components of 
their effective dielectric constants, as will be described. 



D. Debye Random Medium 



The two-point correlation function for a Debye random medium is given by 

^ 2) W-02 = 02(l-0 2 )exp(-r/ 7 ), 



(44) 
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with 7 > 0. This correlation function is plotted in Fig. [2j In the following analysis, we take 
7 = a/2. This correlation function, first proposed by Debye 3l|, was imagined to correspond 
to porous media with cavities of random shapes and sizes (a realization of this medium is 
depicted in Fig. [0(D)). However, it was not known until recently that such a correlation 
function was indeed realizable by a two-phase medium. Using a "construction" procedure, 
Yeong and Torquato 32j demonstrated that the correlation function defined in Eq. (1441) is 
realizable and dubbed such systems Debye random media (see also Refs. Q and jsOj for fur- 
ther discussion of the realizability of Debye random media) . Note that if the two component 
dielectrics and their respective volume fractions are interchanged, the correlation function 
given in Eq. (I44p remains invariant. This property is called phase-inversion symmetry. A 
composite is phase- inversion symmetric if the morphology of phase 1 at volume fraction </>! is 
statistically identical to that of phase 2 at the same volume fraction, 0i 19|, |32j. None of the 



sphere and spheroid dispersions possess this property. It is also important to note that the 
percolation behavior of Debye random media has not been investigated to date. However, 
based on our knowledge of the percolation threshold of the fully-penetrable-sphere model 
and visual inspection of the planar cut through a three-dimensional Debye random medium 
shown in Fig. [TJ we expect that the percolation threshold of the latter is substantially below 
2 = 0.5. 



E. Random Checkerboard 



The random checkerboard model is produced by partitioning space into cubes such that 
each cube is assigned phase 1 or phase 2 randomly according to the volume fraction (see Fig. 
Q^E)). The cubes have side length D = la. The calculation of the radially averaged two- 
point correlation function for this model is not presented here. A detailed derivation of this 
result may be found in Ref. bjj The radially averaged two-point function, S^\r), is plotted 
in Fig. [2j This model may be thought of as being more "ordered' than, for example, the 
Debye random medium because the phases are confined to lie on a grid. In our calculations, 
presented in the next section, we find that the random checkerboard has a greater imaginary 
component, however, than that of hard spheres at volume fraction <p2 = 0.5. This makes 
intuitive sense: in this regime, the hard-sphere system is near crystallization; there is no 
such order-disorder transition for the random checkerboard model. As in the Debye random 
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medium, this model possesses phase-inversion symmetry. For site percolation on a simple 
cubic lattice with nearest-neighbor connectivity criteria, phase 2 would percolate at a volume 
fraction of 02 — 0.312. However, for the related dielectric and conductivity problems, 
it is known that edges and corner points will contribute to the effective properties and 
therefore nearest, next-nearest, and next-next-nearest connectivity criteria must be used. 
This will result in a substantially lower percolation threshold than 02 = 0.312, as one can 
ascertain from the analogous two-dimensional percolation problem on a square lattice. In 
summary, the percolation threshold for the three-dimensional random checkerboard should 
be considerably lower than that for fully penetrable spheres (02 = 0.2895). 



F. Power-Law-Correlated Materials 



The two-point correlation function for these composites are given by 

= (45) 

where the exponent n must obey the inequality n > 3 for the Fourier transform of the 
left side of ( 1451) to exist. The two-point correlation function S\ is plotted in Fig. (j2J) 
with n — 4, which is the value of the exponent that will be considered throughout the rest 
of the paper. A single realization of a power-law-correlated material with this correlation 
function is depicted in Fig. (HF). The degree of correlation in this model is of course 
highly dependent upon the exponent n. Note that as in the cases of the Debye random 
medium and the random checkerboard, the power-law-correlated material is phase-inversion 
symmetric. This correlation function is being put forward in this study for the first time. 
However, it is not necessarily true that any proposed functional form of S2 is physically 
realizable. Any correlation function that corresponds to a realizable material must satisfy a 
number of necessary conditions 30]. Namely, (1) < S% (r) < 02, (2) the radial derivative 
of S\ must be strictly negative at the origin, (3) it must obey the triangle inequality 
S^ir) > S^\s) + S^\t) —02, where r = t — s, and (4) the Fourier transform of S^\r) —cf)\ 
must be everywhere non-negative. We have tested each of these conditions for the power-law 
correlation function, and they are all satisfied. For all other models considered in this work 



have either been known to be realizable or have been recently shown to be 19j,|32j. We choose 



to study power-law-correlated materials because they are reminiscent of structures that are 
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scale-free. While the correlation function given in Eq. (|45j) is, strictly speaking, not scale-free 
(i.e., a function directly proportional to l/r n ), it does asymptotically approach this behavior 
for r ^> a. A purely scale-free function does not satisfy the known realizability constraints. 
Model microstructures that have scale-free correlation functions exhibit interesting clustering 



and percolation properties 19]. This model's percolation behavior will be studied in greater 



detail in a forthcoming study As in the case of the Debye random medium, percolation 
behavior of power-law-correlated materials has not been investigated to date. However, for 
the same reasons given in Section III.D, we expect that the percolation threshold of the 
latter is substantially below 02 = 0.5. 



IV. RESULTS 



In this section, we present results for the effective dielectric constant as predicted by the 
two-point approximation (I52"j) for each of the model microstructures discussed in the previous 
section, i.e., hard spheres, hard oriented spheroids, fully penetrable spheres, Debye random 
media, random checkerboard, and power-law-correlated materials (with exponent n = 4). 
Henceforth, we take £2 > £1, and assume both s\ and 62 are real. We take the wave number 
for a wave propagating through phase 1 to be 27r/(60a). If phase 1 has dielectric constant 

I, then this corresponds to a propagation frequency of 5 GHz. For reasons given in Section 

II. C, we take q = 1 and p = 2 if phase 2 is below its percolation threshold, and we take p = 1 
and q — 2 if it is above its percolation threshold. In the case of the fully-penetrable-sphere 
model, we also evaluate the three-point approximation fl33l) . 



A. Two-Point Estimates 



1. Isotropic Media 

For the hard-sphere model, we employ the two-point approximation given in Eq. fl32|) 
with q = 1 and p = 2 where the spheres comprise phase 2, for 2 < 0.5 (see discussion in 
Section II. C). This choice of two-point approximation is expected to be accurate because 
the equilibrium hard- s phe re model percolates at "jammed" states, which are substantially 



higher than 2 — 0.5 19j. In the case of the fully-penetrable-sphere model, we plot our 



results for a fixed phase contrast ratio but over the entire volume fraction range because in 
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this case we know its nontrivial percolation threshold (see Section III.C). For this model, we 
expect ( 132]) with p = 2 and q = 1 to be valid for volume fractions well below its percolation 
threshold of 02 = 0.2895, and then with p = 1 and q = 2 to be valid well above it. In the 
intermediate region, which is taken to be 0.2 < 2 < 0.4, we interpolate between these two 
regimes using a spline fit in order to approximately account for the fact that this medium 
is crossing its percolation threshold. For Debye random media, random checkerboard, and 
power-law-correlated materials (in which we take the exponent n — 4), the percolation 
thresholds are not known. Therefore, we limit ourselves to analyzing the effective dielectric 
constant of these models for volume fractions that are assumed to be below these thresholds 
(e.g., 02 — 0.1), and well above them (e.g., 02 = 0.5). In the former case, we use the 
two-point approximation f[3"2"]) with p = 2 and q — 1, and in the latter case with p = 1 and 
q = 2. For all of the five isotropic models, we employ the two-point estimate (132]) with the 
expansion given in Eq. (|29|) through third order in k q , and choose the correlation lengths 
to be roughly equal to one another. Different two-point approximations are employed for 
the various statistically homogeneous and isotropic model microstructures below and above 
their percolation thresholds for reasons given at the end of Section II. C. 

For the disordered models described in Section III, we present in Table [I] the two lowest- 
order coefficients of k q a in the expansion of A^\ as given in Eq. (129]) . Note that the zeroth- 
order term is zero, because A% is zero in the isotropic static problem. As can be seen 
from Eq. (|28]) . there can be no first-order term either. The second-order term is necessarily 
real, and the third-order term necessarily imaginary. We call the j th order coefficient in this 
expansion such that = &f\k q ay. 

In Figs. [3] and HI we show the real and imaginary components of the effective dielectric 
constant, respectively, of five isotropic models described in the previous section, plotted 
against dielectric-contrast ratio at a volume fraction 02 = 0.1. Since each model is below 
or assumed to be below its percolation threshold at this volume fraction, the two-point 
approximation ( 132]) is used in conjunction with (129]) in order to calculate A% . The real and 
imaginary parts of the effective dielectric constant at volume fraction 02 = 0.5 are plotted in 
Figs. [5] and [6] respectively. At this volume fraction, all structures except for the hard-sphere 
model are above their percolation thresholds. As shown in Fig. |4j the imaginary component 
of the effective dielectric constant of the hard-sphere model is significantly smaller than 
that of the fully-penetrable-sphere model. This makes intuitive sense because the constraint 
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HS 


FPS 


DRM 


RC 


PLC 




(2) 


(2) 


(2) 


(2) 


(2) 


(2) 


(2) 


(2) 


(2) 


(2) 


02 = 0.1 


0.0512 


0.0304i 


0.0700 


0.0579* 


0.0450 


0.0450* 


0.108 


0.116* 


0.0300 


0.0600* 


02 = 0.2 


0.0658 


0.0279* 


0.120 


0.0990* 


0.0800 


0.0800* 


0.193 


0.206* 


0.0533 


0.1067* 


02 = 0.3 


0.0625 


0.0188* 


0.152 


0.124* 


0.1050 


0.1050* 


0.253 


0.270* 


0.0700 


0.1400* 


02 = 0.4 


0.0512 


0.0107i 


0.165 


0.134* 


0.1200 


0.1200* 


0.289 


0.308* 


0.0800 


0.1600* 


02 = 0.5 


0.0383 


0.0052i 


0.163 


0.130* 


0.1250 


0.1250* 


0.301 


0.321* 


0.0833 


0.1667* 


02 = 0.6 






0.146 


0.115* 


0.1200 


0.1200* 


0.289 


0.308* 


0.0800 


0.1600* 


02 = 0.7 






0.116 


0.0894* 


0.1050 


0.1050* 


0.253 


0.270* 


0.0700 


0.1400* 


02 = 0.8 






0.0770 


0.0575* 


0.0800 


0.0800* 


0.193 


0.206* 


0.0533 


0.1067* 


02 = 0.9 






0.0339 


0.0239* 


0.0450 


0.0450* 


0.108 


0.116* 


0.0300 


0.0600* 



TABLE I: Coefficients of the (k q a) 2 and (k g a) 3 terms in the expansion of (given in Eq. ([29]) ). 
for each of the isotropic models considered here. These are the hard-sphere (HS) model, fully- 
penetrable-sphere (FPS) system, Debye random medium (DRM), random checkerboard (RC), and 
a power-law-correlated system (PLC). For hard spheres and fully penetrable spheres, the spheres 
comprise phase 2. The parameter at? is the j'^-order coefficient. As demonstrated in Eq. (128j) . 
there are no zeroth or first order terms in this expansion. 

that hard spheres may not overlap results in larger spatial correlations (smaller coarseness), 
leading to a smaller imaginary component. Upon increasing volume fraction, the hard-sphere 
model becomes less and less coarse. Close to its freezing point, at 02 ~ 0.5, the resulting 
imaginary component of e e is extremely small, as shown in Fig. [6j At volume fraction 
02 = 0.5, the real part of e e for the hard sphere model is significantly smaller than that 
of the other models, since it is the only one that is not percolating. As Fig. [6] shows, the 
imaginary component of the effective dielectric constant is substantially more sensitive to 
microstructure than the corresponding real component of e e . 
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FIG. 3: (Color online) Real part of the effective dielectric constant of the various isotropic models 
studied here as a function of dielectric-contrast ratio £2/^1 at volume fraction 02 = 0.1 and wave 
number k\ = 27r/(60a). For all models, the two-point approximation f)32|) with p = 2 and q = 1 is 
used since they are each below or assumed to be below their percolation thresholds. In order to 
calculate for each microstructure, Eq. ([29]) is used. 
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FIG. 4: (Color online) Corresponding imaginary part of the effective dielectric constant for the 
isotropic models depicted in Fig. [3j 
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hard- sphere 
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Debye random medium 
random checkerboard 
power- law -correlated 




FIG. 5: (Color online) Real part of the effective dielectric constant of the various isotropic models 
studied here as a function of dielectric-contrast ratio £2/^1 at volume fraction fa = 0.5 and wave 
number k\ = 27r/(60a). For all of the models besides hard spheres, the two-point approximation 
(|32p with p = 1 and q = 2 is used, since they are each above their percolation thresholds. For the 
hard-sphere model, which is below its percolation threshold, we use (132 p with p = 2 and q = 1. 
Eq. ([291) is used to calculate A$ , 
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FIG. 6: (Color online) Corresponding imaginary part of the effective dielectric constant of the 
isotropic models depicted in Fig. [3 On the scale of this figures, the hard-sphere curve is almost 
indistinguishable from the horizontal axis. 
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In Figs. [7] and [8], we plot the real and imaginary components of e e for the fully-penetrable- 
sphere model as a function of volume fraction. The two-point approximation fl32|) is used, 
with p = 2 and q — 1 for < 02 < 0.2 (which is below the percolation threshold of 
02 — 0.2895), and p = 1 and q = 2 for 0.4 < 2 < 1- A spline fit is used to interpolate 
between these two curves to yield an approximation for 0.2 < 2 < 0.4. Eq. f)29p is used to 



calculate A 



(p) 
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fully penetrable spheres 
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FIG. 7: (Color online) Real part of the effective dielectric constant of the fully-penetrable-sphere 
model as a function of volume fraction <p2 of phase 2 at a dielectric contrast ratio £2/^1 = 5 and 
at wave number k\ = 2-7r/(60a). The two-point approximation (|32|) is used with p = 2 and q = 1 
below the percolation threshold of 02 = 0.2895, and p = 1 and q = 2 above it. A spline fit is used 
to interpolate between these two formulas and Eq. (|29p is used to calculate A^f 1 . 
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FIG. 8: (Color online) Corresponding imaginary part of the effective dielectric constant for the 
fully-penetrable-sphere model depicted in Fig. [71 A spline fit is used to interpolate between the 
two formulas. 



Among the isotropic model microstructures that are discussed here, it is possible to carry 
out the full A2 integral as given in Eq. (I28I) . rather than the expansion given in Eq. (I29I) . 
That said, since for both the real and the imaginary parts of , the next term in the 
expansion is smaller by a factor of (k q a) 2 , Eq. ( l29l) gives an extremely good approximation 
to the effective dielectric constant at long wavelength. 

2. Anisotropic Media 

Drawing upon results discussed in Appendix C in which we apply our formalism at the 
two-point level to calculate the effective dielectric tensor of statistically anisotropic with 
azimuthal symmetry, we present results for the axial and planar dielectric constants for an 
equilibrium dispersion of hard oriented spheroids in a matrix for a number of spheroid aspect 
ratios, as predicted by the two-point approximation, given in Eqs. fICll) and flC6j) . Figures 
IH] and [TU] show real and imaginary component results, respectively, for the anisotropic hard- 
spheroid model, giving axial and in-plane dielectric constants of the composite in each plot, 
as predicted by the two-point approximation, given in Eqs. (ICljl and (1C6j) . The former 
plot, Fig. [9] shows that upon increasing the aspect ratio b/a, the axial component of the 
dielectric tensor increases. This makes intuitive sense because in the purely static case, 
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as the microstructure approaches the long-needle limit, the tensor component in the axial 
direction should approach the arithmetic mean (<7i0i + ^fa), which is a rigorous upper 
bound on the effective dielectric constant in the purely static case. We see that the in- 
plane component decreases with increasing aspect ratio. We may again understand this in 
the context of the static regime: in that scenario, the two-point tensor Ag must remain 
traceless. The decreasing in-plane component is a direct result of this property and the fact 
that the axial component increases with aspect ratio. By contrast, we see from Fig. [TOlthat 
both the axial and in-plane imaginary parts of the dielectric tensor increase with increasing 
aspect ratio. 



axial/in-plane component, b/a = 1 

■ ■ ■ ■ axial component, b/a = 2 
— — axial component, b/a = 3 




e 2 /£ l 

FIG. 9: (Color online) Real part of the axial and transverse effective dielectric constants of an 
equilibrium dispersion of hard oriented prolate spheroids, plotted as a function of dielectric-contrast 
ratio £2/^1, for a number of different aspect ratios b/a, as predicted by the two-point approximation, 
given in Eqs. (|C1[) and (|C6|) . The contrast is £2/^1 = 5, the volume fraction is fa = 0.5, and 
k\ = 27r/(60a), where a is the semi-minor axis of the spheroid. 
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FIG. 10: (Color online) The corresponding imaginary part of the axial and transverse effective 
dielectric constants for the same model and cases depicted in Fig. [9l 



Although we do not show results for the oblate case, we have verified that as the aspect 
ratio of the spheroids is increased from the oblate regime (b/a < 1) through the prolate 
regime (b/a > 1), the real part of the axial dielectric constant increases and the real part of 



36. The 



the in-plane dielectric constant decreases, which agrees the behavior given in Ref. 
imaginary components of both the axial and in-plane dielectric constants both increase with 
increasing aspect ratio. 



B. Three-Point Estimates 



Here we apply the isotropic three-point approximation, given explicitly in Eq. (133|) . (with 
p = 1 and q — 2) to the fully-penetrable-sphere model in order to ascertain the importance 
of three-point information. For this model, Table [Til gives the coefficients of (k q a) for the 
first two terms of A% , as given in Eq. ( 13T1) . for a number of volume fractions. We call the 
jth orc [ er coe ffi c i en t in this expansion af \ such that = ^2JL <yf\k q ay . 

The importance of three-point information at high contrast and volume fraction is demon- 
strated in Figs. [11] and [121 which give the (k\a) 2 and (kia) 3 coefficients of the effective 
dielectric constant of the fully-penetrable-sphere model at 02 = 0.5, which is well above 
its percolation threshold. The effective dielectric constant is calculated via the two-point 
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p = 2 


,9 = 1 


p = 1 


q = 2 


(3) 
«0 


(3) 


(3) 
«0 


(3) 


4>„ = o.i 


0.010 


0.012 


0.079 


0.013 


4>„ = 0.2 

' IT 


0.035 


0.040 


0.17 


0.043 


4> p = 0.3 


0.070 


0.075 


0.24 


0.085 


<j) p = 0.4 


0.11 


0.11 


0.31 


0.13 


<j) v = 0.5 


0.14 


0.13 


0.35 


0.17 


<b v = 0.6 


0.17 


0.15 


0.37 


0.20 


4> v = 0.7 


0.17 


0.14 


0.35 


0.21 


4> p = 0.8 


0.15 


0.11 


0.28 


0.18 


4> p = 0.9 


0.10 


0.053 


0.17 


0.12 



TABLE II: Coefficients of the (k q a)° and (k q a) 2 terms in the expansion of (given in Eq. (|3ip ). 
for the fully-penetrable-sphere model, where the spheres comprise phase 2. The parameter a - is 
the j th -ovdev coefficient. As demonstrated in Eq. (|31|) . there are no first- or third- order terms in 
this expansion. 

and three-point approximations, given in Eqs. ( 1321) and ( 133|) . respectively. These equations 
employ the expansions given in Eqs. fl29l) and fl3T|) . For this model, we have taken q = 2 
and p = 1. Clearly, for both the (kia) 2 coefficient (which is real), and the (kia) 3 (which 
is imaginary), three body information plays a significant role at high volume fraction and 
contrast. We also plot the imaginary component of the effective dielectric constant for this 
model, both with and without the third-order contribution in Fig. [121 We see that up to 
a dielectric-contrast ratio of roughly 5, the two-point and three-point approximations are 
in relatively good agreement, but afterwards they increasingly diverge. This serves as a 
test of the convergence of the series. We thus see that up to relatively high contrast ratio 
(i.e., £2/^1 ~ 5), the two-point approximation provides a good estimate, implying that the 
remaining terms in the full series expansion ( fT9l) are negligibly small. 
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FIG. 11: (Color online) The coefficient of (k±a) 2 of the effective dielectric constant as a function of 
dielectric-contrast ratio in the fully-penetrable-sphere model at volume fraction = 0.5, as pre- 
dicted by the two-point and three-point approximations, given in Eqs. ([32]) and ([33]) . respectively. 
We take p = 1 and q = 2 since the system is above its percolation threshold at this volume fraction. 
These equations employ the expansions given in Eqs. (|29p and (|3ip . This term is necessarily real. 
The solid and dashed curves show this quantity with and without the three-point term included 
in the expansion given in Eq. (|17p . 
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FIG. 12: (Color online) As in Fig. ITT1 except for the coefficient of (kia) 
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0.03 



■ — three-point approximation 




FIG. 13: (Color online) The imaginary part of the effective dielectric constant of of the fully- 
penetrable-sphere model at volume fraction fa = 0.5 and wave number k\ = 27r/(60a), as predicted 
by the two-point and three-point approximations, given in Eqs. (|32p and (|33p . respectively. We 
take p = 1 and q = 2 since the system is above its percolation threshold at this volume fraction. 
These equations employ the expansions given in Eqs. (l29|) and (j3~Tj) . The solid and dashed curves 
show this quantity with and without the three-point term included in the expansion given in Eq. 

CZD. 



V. CONCLUSIONS 



We have derived exact strong-contrast expansions (1191) for the effective dielectric tensor 
e e of electromagnetic waves propagating in a two-phase composite random medium with 
complex-valued isotropic components in the long-wavelength regime. These expansions are 
not formal but rather are explicitly given in terms of certain integrals over the n-point 
correlation functions that statistically characterize the medium. To our knowledge, such an 
exact representation has not been given explicitly before. The nature of the strong-contrast 
expansion parameter results in a radius of convergence of the series (1191) that is significantly 
widened beyond that of a weak-contrast expansion (i.e., the simple difference of the dielectric 
constants of the two phases). Because the expansions can be considered to be perturbations 
about the solutions of the dielectric tensor of certain optimal structures, we argued that 
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the first few terms of the expansion ([19]) should yield a reasonable approximation of e e , 
depending on whether the high-dielectric phase is below or above its percolation threshold. 
In particular, truncations of the exact expansion led to two- and three-point approximations 
for e e , which we applied to a variety of different three-dimensional model microstructures, 
including dispersions of hard spheres, hard oriented spheroids and fully penetrable spheres 
as well as Debye random media, random checkerboard, and power-law-correlated materials. 

We payed special attention in our analysis to case in which the components have real di- 
electric constants but where the effective dielectric tensor possesses imaginary components 
due to disorder, a phenomenon which is essential to applications such as remote sensing (e.g. 
in the calculation of the backscattering coefficient ll|. In examining two-point approxima- 
tions of the effective dielectric constant for statistically homogeneous and isotropic media 
with component phases having purely real dielectric constants, we found that the imaginary 
part of e e is related to the coarseness for very large windows, which may be regarded to be a 
crude measure of disorder. Among other results, we found that the equilibrium hard-sphere 
model for volume fractions up to its freezing point exhibit a much smaller imaginary com- 
ponent of the effective dielectric constant than the other four statistically homogeneous and 
isotropic model microstructures studied here. 

For dispersions of fully penetrable spheres, we analyzed the behavior of the effective 
dielectric constant using the two- and three-point approximations, Eqs. fl32l) and fl33l) . 
respectively. Our results suggest that truncation of the exact expansion (fl9l) at the two-point 
level yields good convergence up to relatively high phase-contrast ratios (~ 5). However, 
the two approximations increasingly diverge from one another for higher values, showing the 
importance of using three-point information at sufficiently high contrast ratios. 

In order to demonstrate the application of our formalism to a statistically anisotropic 
and hence macroscopically anisotropic media, we examined dispersions of equilibrium hard 
spheroids. It was shown that as the aspect ratio b/a was increased, the real parts of the axial 
dielectric constant and the in-plane dielectric constants increased and decreased, respectively. 
The imaginary components of the axial and in-plane dielectric constants were both found 
to increase upon increasing the aspect ratio b/a from 1 (i.e., sphere point). 

The dichotomy between periodic media, which have zero imaginary component in their 
effective dielectric constants, and disordered media, which do have such a component, begs an 
important question: which is true of quasiperiodic structures? These structures possess long- 
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range order but have no translational symmetry, are therefore are in a sense intermediate 
between the latter two regimes. In Appendix B, we show that at the two-point level, an 
imaginary component is obtained only if the Fourier transform of the correlation function 
is nonzero at k q . For two-phase media with periodic structure (i.e., crystals), this implies 
that at this level there is no imaginary component, but this argument does not hold for 
media with quasiperiodic structure (i.e., quasicrystals) |40|. Since the diffraction pattern for 
quasicrystals possess only discrete Bragg peaks, the coarseness is necessarily zero. However, 
this does not imply that the imaginary component of e e is identically zero, since higher-order 
k q coefficients may contribute. The calculation of the effective properties of quasicrystal two- 
phase media and a fundamental understanding of their wave propagation properties remains 
a challenging open question. 

In future work, we plan on applying the strong-contrast formalism of this paper to a 
number of new model microstructures and explore other comparison materials in the spirit 



of Ref. 



271 to yield even better approximations for e e . The procedure presented here is 



applicable in any dimension and is thus well suited to the study of scalar Helmholtz equation 
for arbitrary space dimension. One application of this case lies in the calculation of effective 
properties of phononic systems. Another possible extension involves generalizing analogous 



elastostatic results 



19 



391 ] to the dynamic case. 
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APPENDIX A: TWO-POINT EXPANSION FOR ARBITRARY COMPARISON 
MATERIAL 

The strong-contrast expansion ffl9|) for the effective dielectric tensor was derived for the 
choice of the comparison material such that Eq = e q (q = 1 or 2). This this choice sig- 
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nificantly simplifies the analysis and has the desirable feature that it can be regarded to 
be a solution that perturbs around the effective dielectric tensor for certain optimal mi- 
crostructures. However, for other microstructures, different comparison materials may offer 



advantages in terms of better series convergence |26| and more accurate approximations 



27| . Here we present an equivalent relation to Eq. (1321) for macroscopically isotropic media 



but for an arbitrary comparison material with dielectric constant Eq: 



PeD P Q + (PpO ~ P q0 )<Pp (f3 g0 + ((3 p0 ~ P q o)<P P f 

where /3 eQ _ i = (e e - e )/(e e + 2e ), p p0 = {e p - e )/{e p + 2e ), and (3 q0 = (e q - s )/(£ q + 2e ). 

£o is taken to be the Bruggeman effective dielectric constant, but this is not the 
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In Ref. _ 

best choice for general microstructures [ill]. For example, for the large class of dispersions 
discussed in Section II. C, the choice e = e q is better. 



APPENDIX B: PROOF THAT A 2 IS PURELY REAL AT THE TWO-POINT 
LEVEL 

All of the random media studied here have had dielectric tensors with non-zero imag- 
inary component. Thus, electromagnetic waves propagating through these materials will 
attenuate, albeit with small decay constants. Roughly speaking, the physical cause for 
this attenuation is the fact that the waves are scattered off of the heterogeneities, and the 
scattered waves are no longer coherent with the propagating wave, and thus, this energy 
is "lost", when homogenization is applied in Eq. ffT^l) . However, this is not so for peri- 
odic heterogeneous media. In the language of solid-state physics, scattering does not take 
place because waves are allowed to propagate coherently as Bloch waves (as opposed to in 
the form of pure plane waves). In this Appendix, we show this explicitly at the two-point 
level. Although a single periodic configuration is statistically inhomogeneous, our formu- 
lation (valid for statistically homogeneous media) can still be applied by first performing 
a translational average, i.e., averaging over uniformly random displacements of the origin. 
This produces averaged quantities, such as the correlation function defined by Eq. fl25l) . 
translationally invariant. Thus, invoking an ergodic hypothesis [lj]], the ensemble average is 
equal to an infinite-volume average over the variable ri of a single (periodic) realization. We 
limit ourselves here to periodic media that yield macroscopically isotropic dielectric tensors. 
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The two-point coefficient (1251) in the macroscopically isotropic case can be written as 
follows: 



(4tt) J r 



(Bl) 



We define a "difference indicator function" as follows: 

V(r) -<j> p (B2) 

such that (^ p) (r) - 0j) = (V(R)V(R + r)), where the average is taken over the dummy 
variable R. 

Representing this real-space correlation function in its spectral form, we have 

(V(R)V(R + r)> = J j^-V(q)V(-q) exp(*q • r). (B3) 

where V(q) is the Fourier transform of the difference indicator function at wave vector q. 
Inserting this expression into Eq. (1B1I) . we have 

.00 _ 2k l f ,y kqr f d( i i> 



^ = (4$ 7 y (^(^(-qWCiq-r) 

2k 2 f da ~ ~ f e ik " r 

= W)j (2^ V(q)V( - q) / ^ ex PN-)^- (B4) 

We see that the inner integral is just the Fourier transform of the Green's function of the 
Helmholtz equation. This inner integral is thus simply Air/(q 2 — kT). Thus, taken together, 



we have 



Af = ^ / ^V^H-d)^ (B5) 



(4tt) J (2vr) 3 v q 2 -k 



We can see that there are two poles in the last integral located at q = ±fc 9 . Clearly, for 
a periodic medium, V r (q)V r (— q) is identically zero at every point, save the reciprocal lattice 
vectors, where there are delta functions. If we assume long wavelength, then k q G, where 
G is the magnitude of the smallest non-zero reciprocal lattice vector. Thus, the pole will 
not enter in to the calculation! This proves that there can be no imaginary component to 
this integral, and this must be the case for all periodic structures. For the case of a random 
medium, since there are indeed correlations at long wavelength (of which there are of course 
none in the periodic case), there is a non-zero residue at the pole, and the integral can be 
complex. 
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APPENDIX C: FOR AZIMUTHAL SYMMETRY 



We now discuss how to explicitly express the tensor two-point coefficient [cf. (I2"T1) ] 
for statistically anisotropic media with an axis of symmetry (say, the z axis) i.e., azimuthal 
symmetry. An example of such a microstructure is a dispersion of hard oriented spheroids, 
which we discussed in Section III.C. Here we follow the methodology of Torquato and Lado 



36( | . who evaluated similar integrals for the purely static case. If we align the Cartesian 
coordinate system with the principal axes frame, the two-point coefficient A^ is diagonal, 



i.e., 



( U ^ 



A 



(p) 



V 



(Cl) 



/ 



U 
V 

where U is the in-plane (i.e., x-y plane) component and V is the axial (z axis) component. 
We may explicitly write the two-point tensor as 



4^ = — I drexp(ik q r 



(3 



3ik q r 



k 2 r 2 )rr 



-1 + ik q r + k 2 r 2 )I 



S. 



ip), 



(C2) 



where we have inserted the dyadic Green's function explicitly. Retaining terms through 
order k q , gives the following explicit expressions for U and V in spherical coordinates: 



U 



3 

Arc 



*(-l + £sm 2 (<S9) 



and 
V 



2ik\ 3 



dr 



dr 



S. 



ip). 



s. 



(p), 



+ 



k q 3 



^(1 + Isin 2 ^)) 



S. 



(p). 



(C3) 



47T 



f- / =!(-! + 3cos a (e))(5 a w (r 



2 4tt 



-(l + cos 2 (£)) 
r 



r) 



2?fc| 3 



dr 



ip), 



(C4) 



3 4vr 

where we have used the fact that the two in-plane matrix components must be equal in 
order to simp 



In Ref. 
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ify the integrals. 

the first integral in the expression for V was evaluated analytically for the 
case of hard oriented spheroids. In what follows, we specialize to this model. We evaluate 
the first two integrals in each of the relations (IC3I) and flC4j) numerically. Note that the 
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third integral in each expression, which we denote by U3 and V3, respectively, are identical, 
and may be obtained analytically as follows: 



U« = V* 



2ik\ 3 

Anile* 3 
3 4^ 



dr 

1 

1 



S%(r;b/a)- 



1 w r 

dcos(9) / drr 2 S$% s {<T (r/<T{9)); 1) - <j>; 



ik* 



dcos(9) 



ikl / dcos(9) 



-1 



a 



(0) 





(irr 2 
3 



S^ s {<r/a{9))- 1 l)-<t> : 



drr' 



2iki 



drr' 



(C5) 



where we have used the mapping between the two-point function S^ HS (r; b/a) for the hard- 
spheroid model and S2 #s(r; 1) for the corresponding hard-sphere system described in Section 
III.C. Note that the radial integral in the last line of (1C5|) is proportional to the zero-wave- 
vector of the structure factor of the hard-sphere system. For an equilibrium distribution, this 
may be obtained analytically directly from the Percus-Yevick approximation, for example. 
The two-point approximation to the dielectric tensor is obtained by truncating the series 
expansion given in Eq. ( 1191) after the second term to yield 



P 2 pq <Pl(e e -e q I)- 1 (e e + 2e q I) 



(C6) 
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